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ABSTRACT 

As the Galaxy evolves, the abundance of deuterium in the interstellar medium 
(ISM) decreases from its primordial value: deuterium is "astrated". The deu- 
terium astration factor, /d, the ratio of the primordial D abundance (the D to 
H ratio by number) to the ISM D abundance, is determined by the competition 
between stellar destruction and infall, providing a constraint on models of the 
chemical evolution of the Galaxy. Although conventional wisdom suggests that 
the local ISM (i.e., within ~ 1 — 2 kpc of the Sun) should be well mixed and 
homogenized on timescales short compared to the chemical evolution timescale, 
the data reveal gas phase variations in the deuterium, iron, and other metal 
abundances as large as factors of ~ 4 — 5 or more, complicating the estimate 
of the "true" ISM D abundance and of the deuterium astration factor. Here, 
assuming that the variations in the observationally inferred ISM D abundances 
result entirely from the depletion of D onto dust, rather than from unmixed ac- 
cretion of nearly primordial material, a model- independent, Bayesian approach 
is used to determine the undepleted abundance of deuterium in the ISM (or, 
a lower limit to it). We find the best estimate for the undepleted, ISM deu- 
terium abundance to be (D/H)ism > (2-0 ±0.1) x 10~ 5 . This result is used to 
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provide an estimate of (or, an upper bound to) the deuterium astration factor, 
/ D = (D/H) P /(D/H) ISM <1.4±0.1. 



Introduction 



Deuterium is create d in an astrophysically intere s ting; abundance only during big bang 
nucleosynthesis (BBN) ( IBoesgaard fc Steigmanl Il985l ; ISteigmanl 120071 ) . after which, in the 
post-BBN Universe, its abundance (?/d = 10 5 (D/H)) decreases monotonically due to the 
processi ng of gas through su c ceeding generations of star s where deuterium is completely de- 
stroyed (jEpstein et al.lll976l ; IProdanovic fc Fields! 120031 ) . Consequently, deuteri um plays a 



special role in cosmology, nuc 
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Prodanovic fc Fields! 2008 ). Its relatively simple Galactic evolution permits us to use 
deute rium to determine the fracti on of interstellar gas that has been processed through 
stars (ISteigman fc Tosil Il992l Il995l ). By comparing the primordial and Galactic deuterium 



abundances on e can learn about and discriminate among different Galact i c chemical evo- 



lution models flVangioni-Flam et al.l Il994t iTosil Il996t iRomano et al.l l2006t ISteigman et al. 



20071 : IProdanovic fc Fieldsl I2OO8 L Together with non-BBN constraints on the baryon (nu 
cleon) density from observations of the cosmic microwave background ( jSpergel et al.l 12007 



Dunkley et al.l 12009 : Komatsueta 



is predicted by BBN dCvburt et al 



2009 



2010) , the primordial deuterium abundance, t/dp 



200.4 ICoc etldlhooi ISteiginan&Tl: Isimha fc Steigman 



2008 ; Cyburt et al. 20081 ). The predicted abundance is in excellent agreement with observa- 
tions of deuterium in high-redshift , low-m etallicity QSO absorption line systems (QSOALS) 
Jc-'Meara et al.ll2006[ Fettini et al.ll2008l ). 



Since deuterium is destroyed in the Galaxy as gas is cycled through stars, (D/H) ISM < 
(D/H)p. However, all successful chemical evolution mod els require infall to the disk of 
the G alaxy of unprocessed (or, nearly unprocessed) gas (ITosil 1 19961 ; IProdanovic fc Fields 
20081 ). and such deuterium-rich (and metal-poor) gas would raise the ISM D/H ratio closer 
to the primordial value. As a result, comparison of the primordial and ISM deuterium 
abundances provides a constraint on infall and, therefore, on chemical evolution models. 
Observations over the past decade and more of the deuterium abundance in the relatively 



local interstellar medium (ISM) reveal an unexpectedly large scatter in D/H (IJenkins et al. 



19991 ; ISonneborn et al.l |2000| ; iHebrard et al.l |2002| ; iHoopes et al.l 120031 ) , challenging the con- 
ventional wisdom of a well mixed ISM. For example, as shown in Figure HJ the absorption 
line measurements from the Far Ultraviolet Spectroscopic Explorer (FUSE) reveal variations 
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Fig. 1. — The logs of the deuterium abundances versus the logs of the Hi column densities 
[cm -2 ] for the 49 FUSE LOS (see the text). The filled symbols are for the 41 LOS which have 
iron abundance data, while the open symbols are for the 8 LOS which lack iron abundances. 
The squares (blue) are for the LOS within the Local Bubble (LB) and the circles (red) are 
for the non-LB (nLB) LOS. The solid line is at the mean D abundance for the LB LOS (log 
2/d,lb = 0.19); the dashed line is its extension to the nLB LOS. 



of a factor of ~ 4 (0.5 < ?/d,ism = 10 5 (D/H) ISM ^ 2.2) in the gas-phase D/H ratios over 
lines of sight (LOS) to background stars within ~ 1 — 2 kpc of the Sun. Moreover, the 
variations in the observed, gas phase D/H abundanc es are found to correla t e positively with 



the a bundances of refractory elements such as Ti (IProchaska et al.l 120051 ; lLallement et al. 



200 8 ). and a similar c orrelation is also fou nd between D and other metals such as Fe and 
O (ILinsky et al.ll2006l ; ISteigman et al.ll2007l ). However, this positive correlation between the 
abundances of D and the metals is opposite to the trend expected from stellar nucleosyn- 
thesis (D decreasing as the metals increase). Motivated by the observed correlations and by 
the very large spread in the D/H ratios inferred from the FUSE and other data sets, it has 
been proposed that the large variat ions in gas-phase D/H may be d ue to the depletion of gas 
phase deuterium onto dust grains (1 Jural Il982l ; iDraind 12004 120061 ). If this is the case, then 
the FUSE (and other) gas-phase absorption-line measurements reveal that depletion has not 
been well mixed (homogenized) in the ISM and, the data may only provide a lower limit to 
the true (undepleted) ISM deuterium abundance and, therefore, only an upper limit to the 
deuterium astration factor, / D = ?/dp/2/d,ism- 

In Figure [TJthe logs of the deuterium abundances (log y D = 5 + logN(Di) — logN(Hi)) 
are shown as a function of the logs of the H i column densities for the 49 LOS with D i column 
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densities from Table 2 of 



( 120061 ) and iDupuis et al 



inskv et al.l ( 120061 ) . supplemented by data from lOliveira &: Hebrard 



( 120091 ). Of these 49 LOS, 41 have iron abun dance measur e ments 



(filled symbols); 21 of the 49 LOS are within the Local Bubble (LB; see iLinsky et al.l (120061 ) 
for a discussion of the LB); the remaining 28 non-Local Bubble (nLB) LOS are toward stars 
beyond the LB (see §2.1). 

The unexpectedly large spread among the observational ly inferred ISM D a bundances 
complicates any estimate of ?/d,ism- Recognizing this point, ILinsky et al.l ( 120061 ) chose for 
their estimate of a lower bound to the true (undepleted) ISM deuterium abundance the mean 
of the five highest D/H ratios finding ?/d,ism > 2.17 ± 0.17 or, when including corrections 
(see the discussion in §2.1) for N(Hi) and N(Di) for those lines of sight outside of the Local 
Bubble, ?/d,ism > 2.31 ± 0.24. These estimates are quite close to th e lower bound to t he 
primordial abundance estimated from the QSOALS, ?/dp = 2.821^ ( Pettini et al. 2008l )R 



suggesting a small upper bound to the D astration factor, /p < I.SO^q^ or, an even smaller 



value, /d < 1-22 ± 0.15 for the LB-corrected, nLB deuterium abundances. To account for 
such a high ISM deuterium abundance and such a small D astration factor, a very high infall 
rate of pristine material would be needed, challenging many Galactic chemical evolution 
models (jllomano et al. 20061; Prodanoyic &: Fields 2008 ). By limiting themselves to the five 
highest D abundances, ILinsky et al. (bood ) ignore the lower deuterium abundances along 
many more LOS which are consistent with them within the errors, potentially biasing their 
estimates of 2/dism and of /p- 



In an attempt to address this issue, ISteigman et al.l (120071) (SRT) , used the 18 highest 
D/H ratios from the FUSE data (see Table 3 of ILinsky et al.l (120061 )). finding an ISM D 



abundance of ?/d,ism = 1-88 ± 0.11, corresponding to a D-astration factor /b < I-SO^q;}^, 
consistent with at least some of the otherwise successful chemical evolution models discussed 



in SRT. In fact, the data in Table 2 of ILinsky et al.l (120061 ) reveal that there are 19 LOS 



1 Since this estimate of j/dp relies on only 7 high redshift, low metallicity LOS and, since the dispersion 
in deuterium abundances among them is unexpectedly large, some prefer to adopt a so-called "WMAP D 
abundance" . WMAP does not observe deuterium. Rather, the WMA P-determined ba r yon d ensity parameter 
may be used in a BBN code to predict the relic D abundance. If the iKomatsu et al.l (|2010h estimate of the 
baryon density is adopted, the BBN predicted primordial D abundance is z/dp = 2. 5 ± 0.1. However, the 
WMAP collaboration also provides an estimate of the effective number of neutrinos, (|Komatsu et al.ll2010l ) 
which, when used along with their baryon density estimate in a BBN code, leads to a different predicted 
primordial D abundance yop = 3.0 ±0.4. The difference between these two predictions reflects the difference 
between the standard model effective number of neutrinos expected and the WMAP observed value, which is 
within ~ 1.5cr of expectations. So, which, if either, "WMAP D abundance" is preferre d? Here, we compare 



observations to observations, not to mo del- dependent predictions and, we adopt the iPettini et al 
value for yup- 



d2008l) 
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with central values of log yn > 0.20. The weighted mean (along with the error in the mean) 
for these 19 D abundances is log ?/D,i9 = 0.26 ± 0.01, corresponding to yn,i9 = 1-8 ± 0.1. 
For these 19 LOS the reduced x 2 — 0.85, confirms that the weighted mean provides a good 
description of their D abundances. As more LOS with lower D abundances are added, the 
weighted mean decreases, but the reduced x 2 increases, so that ?/d,ism > 1-8 ±0.1 is likely a 
robust lower bound to the ISM D abundance and, /d ^ 1.5 ±0.1 (log /d < 0.19 ± 0.03) is a 
robust upper bound to the deuterium astration factor. 

Surely, there must be a better way to find a reliable estimate of the maximum value of 
the deuterium abundance in the local ISM while accounting for the observational errors in the 
individual D abundance determinations. In §2 we describe a Bayesian analysis designed to 
find the best estimate of the maximum, gas phase (undepleted), ISM deuterium abundance, 
Z/D,max; from data with non-negligible errors and we apply it to the FUSE data set. On 
the assumption that the spread in the observed D abundances is the result of incompletely 
homogenized dust depletion, ?/d,ism > 2/D,max and / D = 2/dpA/d,ism < 2/Dp/2/D,max- Our 
results are summarized and our conclusions presented in §3. As already noted, the dust 
depletion hypothesis suggests that there should be a correlation bet ween the D and F e 
abundances. Although this trend appears to be present at some level (ILinsky et al.ll2006l ). 
a deeper analysis sug gests that the simples t interpretation of the depletion hypothesis may 
need to be modified (ISteigman et al.l 120071 ). For example, since D is likely attached to the 
mantles of dust grains ( Frielens 1983 ). while Fe is mostly contained in the cores of the grains, 
D will be removed from grains more easily than Fe. Thus, although strong shocks may 
destroy grains completely , weak shocks may only remove D from grains while Fe might stay 
locked within their cores ( jDraine fe Salpeterlll979l ). This effect then may explain the scatter 
observed in the relation between the D and Fe depletions. Thus, in a follow-up paper in 
preparation, the consistency of the depletion hypothesis with the FUSE data for D and Fe 
(and O) is investigated, and so me possible modifications to its simplest version are explored 
( ISteigman fe Prodanoviq|2010l ). 



2. A Bayesian Analysis Of The ISM D Abundances 



To avoid imposing any prior prejudice on which LOS should be included and which 
excluded in our analysis or, on how the observed D abundances may, or may not correlate 
with iron (or other metals), we adopt a statistical, model-independent method for determin- 
ing the undepleted Galactic deuterium abu ndance. Our approa ch follows closely the model- 
independent Bayesian analysis developed bv lHogan et al.l (|1997) t o deter mine the primordial 
helium abundance. It is useful to recall the problem iHogan et al.l ( 119971 ) addressed and how 
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they solved it. Helium-4 is produced abundantly during BBN and in the post-BBN Universe 
its abundance is enhanced by stellar produced 4 He. As a result, the 4 He mass fraction, Y P , 
is expected to increase from its primordial value, Yp, along with the metallicity, Z. Extrapo- 
lation of the observed Y versus Z relation to zero metallicity results in an estimate of Y min = 
Yp. The problem for 4 He is that the form of the Yp versus Z relation is a priori unknown 
and, may even differ from object to object (low metallicity, extragalactic Hn regions). This, 
in combination with the errors in the observed val ues of Yp and Z, com plicates the derivation 
of Yp from the data. The challenge con 



ronting IHogan et al.l ( 119971 ) was to identify Y 



The Bayesian approach they developed (IHogan et al.lll997l ) is described below. If, indeed, 
the observed spread in D abundances results from dust depletion, then their 4 He problem is 
entirely analogous to the one we confront in using the D abundance data, with its errors, to 
infer the maximum, gas phase D abundance. 

Here, the FUSE set of Galactic ISM deuterium observations is analyzed assuming only 
that there exists a "true" , uniform, ISM deuterium abundance whose gas phase value may 
have been reduced by the depletion of D onto dust grains, with the amount of depletion 
possibly varying from LOS to LOS. The a priori unknown distribution of depletions is char- 
acterized in terms of a Bayesian probability distribution (a "prior"), and the data themselves 
are used to determine both the true ISM abundance ?/d,ism = Z/D,max or, a lower bound to it, 
2/d.ism > 2/D,max, as well as a measure of the amount of depletion, the depletion parameter 
w = yD,max — 2/D,mm- No prior assumptions about which LOS may have been affected by dust 
depletion are imposed and the entire data set is analyzed in an unbiased, Bayesian manner. 



2.1. The Data 



In our analysis we use the FUSE ISM deuterium abundance data for 46 LOS from 
Linsky et al.l (120061). together with three, more recent measurements tow ards HD41161, 
HD53975 Joiiveira fc Hebrar J2006f ) and REJ1738+665 JDupuis et al.ll2009h. Of t he 49 LOS 



with deuterium abundance measurements, 21 are within the LB (see lLinsky et al.l (120061 ) for 
a discussion of the LB); the remaining 28 are nLB LOS, towards stars beyond the LB. While 
the star 31 Com is more distant than all the LB stars, the LOS to it has a very low H i column 
density and an ave rage Hi volume densi ty (n (Hi)=N(Hi)/d) m uch smaller than for all LB 
LOS. According to Piskunov et al. ( 1997 ) and Dring et~aU ( 1997 ). the absorption feature to- 
ward 31 Com is at a velocity which is inconsistent with the LB and this LOS likely lies within 
the hot, ionized, Local Bubble (which may account for the "high" D and Fe abundances) in 
the di rection of the north galactic pole. For a contr ary point of view, s ee iRedfield fc Linsky 
(120081 ). In contrast to the assignment in Table 2 of iLinsky et al.l (120061 ). we include 31 Com 
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along with the nLB LO S (see Figure |T). Ho wever, this assignment has negligible impact on 
our quantitative results. iLinsky et al.l (120061 ) also list D/H ratios for the nLB LOS corrected 
by them for assumed average foreground LB D i and H i absorption. If velocity information 
were available it could be used to separate foreground absorption in t he LB from that due to 
gas lying beyond the LB. In the absence of such data, ILinsky et al.l (120061 ) assume that the 
LB extends to N(Hi)lb = 10 192 and multiply N(Hi)lb by an adopted average LB D/H ratio 
to find N(D i)lb- These two average column densities are subtracted from the D i and H i col- 
umn densities observed for the individual nLB LOS and the ratio of the corrected column 



densities is used to find corrected nLB D/H rati os. ILinsky et al. (2006) base their assumed 
LB Hi column density on the Nai observations of lLallement et al.l ( 120031 ). who find tentative 
evidence for a LB "wall" of cold dense gas with N(Hi) ~ 10 19 5 . This procedure can bias the 
"corrected" nLB deuterium abundances. For example, it enhances the D abundances along 
those nLB LOS where the observed D abundance exceeds the adopted LB D abundance and 
decreases the D abundances for those LOS where the observed D abundance is less than the 
adopted LB value. The magnitude of the correction increases with the value of the "average" 
Hi column density adopted. In fact, the observed LB Hi column and volume densities are 
distributed very inhomogeneously, varying by nearly a factor of 30 along different LB LOS. 
This results in a scatter plot for the column densities as a function of distance to the back- 
ground star. Even for the subset of the most distant LB stars with al ~ 70 — 80 pc, N(Hi) 
varies by a factor of ~ 16, 10 17 9 < N(Hi) < 10 191 . The structure of the LB is very complex 
(ILallement et al.l 120031 ) and the location of the LB "wall" varies fr om 65-150 pc. Thus, the 
value for the LB Hi column density adopted by ILinsky et al.l (120061 ) may be an overestimate, 
leading to an overestimate for their inferred value of ?/D,max, since a smaller choice for the 
average foreground LB Hi column density will result in a smaller correction to N(Di) and a 
lower estimate of j/D,max- For these reasons, in our analysis, we use the uncorrected Hi and 
Di column densities for all nLB LOS. 



Following ILinsky et al.l (120061 ). we adopt as a working hypothesis that the large scatter 
in the observed, gas phase, ISM deuterium abundances is a reflection of the preferential 
depletion of deuterium (relative to hydrogen) onto dust grains. Therefore, it is assumed 
that in the ISM deuterium has a total LOS column density N(D)ism = N(D) gas + N(D)d us t5 
where the sum includes the observed, gas phase component, N(D) gas = N(Di) bs ; an d an 
unobserved, dust-depleted component, N(D)d U st- Along any LOS, N(D)ism > N(D) gas . The 
data suggest, and we assume, that hydrogen is negligibly depleted onto dust and that the 
fraction of H tied up in molecules in the diffuse ISM probed by FUSE is usually negligible 
(when H 2 is observed, its column density is included in the budget for gas phase hydrogen), 
so that N(H)i SM = N(H) gas . 



On the assumption of a uniform D abundance (gas plus dust) in the local ISM, the 
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dispersion among the observed gas phase D abundances, yr>,i-, reflects the observational 
errors in N(D i) and N(H i), along with any spatial variations in D depletion. Since deuterium 
may be depleted along ALL local ISM LOS, yo,? < 2/D,max, where yv>,max is the maximum 
gas phase (undepleted) deuterium abundance, consistent with the observational errors, and 
2/d,ism > 2/D,max- Within the local ISM there will be a maximum depletion so that yn,% > 
2/D,min- Along the zth LOS, the observed, gas phase D abundance, yr>,i (yB,i = yD,t,obs)> ma y 
differ from the maximum undepleted abundance by an amount Wi where, Wi = ?/D,max — Vd,u 
and < Wi < w = yn,ma,x — yr>,mm- m the absence of observational errors, the observed 
D abundance along the ith LOS would lie between ?/D,min and t/D,max and would reflect 
the difference between the maximum undepleted ISM abundance jto,max> and the spatially 
varying value of w, Wi — yr>,max — yv>,i- Observational errors complicate the task of using the 
data, {yn,i}, to identify y D ,max (and y Djmin or, w). 

Since real data do have errors, the observationally inferred D abundance along the 
ith LOS, yr>,i, will differ from the true gas phase D abundance, yB,i,T, by an amount 5, 
{yD,i = y~D,i,T + Si). For an individual measurement the difference between the observed and 
true D (gas phase) abundances, Si, is unknown. If it is assumed that Si is a random variable 
drawn from a zero- mean Gaussian of width cij, then the probability distribution for y^^, 
given a true gas phase abundance 2/d,«,t along the ith LOS is 

P(yD ti \yn,i,T) = -y^e-^-^T' 2 ^ (1) 

V27rcr i± 

In eq. (II]) we allow for asymmetrical measurement errors so that o"j + corresponds to yrj,i > 
y Dii)T and to y Dji < yr>,i,T- 



2.2. Bayesian Formalism 



For a set of data with errors {yn,i, Si}, the Bayesian approach described by lHogan et al. 



( 119971 ) enables a determination of two parameters, ?/D,max and w (or |/D,min)- Along each LOS 



in the local ISM the true abundance is related to ?/D,max by y^> % t = yv> max _ w 



The likelihood of finding particular values of yo.max and of w (or yD,min)> =^(l/D,max; w), 
from a sample of measurements yv ;i , with errors Si is, 

^(yD,m ax ]w) = Y[P(yD,i\yD,m ax ;w) (2) 

i 

;w) (3) 
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where P(jto,i|$to,i,T) is the probability distribution from eq. (JT]), relating each observed abun- 
dance ?/D,i to its underlying true value yn,i,T- P (2/D,i,T||/D,niaxi w ) is the dust depletion proba- 
bility distribution - the probability of finding the true, dust depleted, LOS deuterium abun- 
dance ?/D,i,T) given the values of the maximum and minimum ISM deuterium abundances, 
2/D,max and y D ,min = 2/D,max - w respectively. Our probability distribution, P(?/D,i,T|z/D,max; w), 
determines the distribution of the Wi. To avoid any a priori model dependence and to limit 



ourselves to the fewest assumptions, following iHogan et al.l (119971 ) we initially adopt the 



simplest form for this probability distribution - a top-hat function: 

D/ I \ J l l W > 2/D,min < VD,i,T < 2/D,max , 

Pfo^Vn^w) = ^ Q ; otherwige (4) 

where yn,mm = Jto.max - w. The probability distribution is normalized to unity when inte- 
grated over yn,i,T- Integrating over the range of possible deuterium abundances for N data 
points, we find the likelihood distribution for the maximum gas phase deuterium abundance 
in the local ISM, yT> max , and for w, 



cs, ( TT ' t ( V~D,i — Jto.miiA t ( VD,i ~ ?/D,max\ 

^o,„„ = erf {~^ r ) - erf {^~) 



(5) 



Equation (jSJ) is evaluated numerically for a range of ?/D,max and w (or yn,mm) values to find 
the combination of these two parameters which maximizes the likelihood. In this way the 
set of deuterium observations, including the errors, is used to find the most likely values of 
2/D,max and the depletion parameter, w, avoiding any prior assumptions about which LOS 
may, or may not, be affected by depletion of deuterium onto dust or, about how the abun- 
dance of D may, or may not, be correlated with the abundance of iron or other metals. 



In addition to adopting the top-hat depletion probability distribution, following IHogan et al 



( 119971 ) we also explore the consequences of choosing two other, asymmetric probability dis- 



tributions: 



r>l I , I %D,i,T - 2/D,min)/w 2 , Z/D,min < 2/D,i,T < Z/D,max , R v 

P{v»«r\Vo^w)={ Q otherwige (6) 



n/ I \ j 2 (|/D,max - 2/D,i,T)/w 2 , 2/D,min < 2/D,i,T < Z/D,max 

P(y^My^w) = | Q ^ otherwige (7) 

where eq. (jUJ) is a positive bias distribution, favoring smaller depletion of deuterium, while 
eq. (J7J) is a negative bias distribution, favoring larger D-depletion. At the Referee's suggestion 
we have also analyzed the effects of another probability distribution - a bimodal, M-shaped 
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bias distribution defined in eq. OH]), which favors both low and high depletion of deuterium, 
while moderate depletion is disfavored. Finally, we also explored the complementary, A- 
shaped probability distribution presented in eq. OH])- Both additional distributions are defined 
with respect to the midpoint y M = (2/D,max + Z/D,min)/2 = ?/D,max - w/2. 



-P(2/D,i,T|Z/D,max; w) 



»(HM+»-KD,max) ' ^D,min < Jto.i/T S 2/M 

2(j/M-yp,i,T) ,, ^ 



^(|/DAT|2/D,max,W) - <; 2(y D . < , T -y D .; nllx ) (9) 

All probability distributions are normalized to unity when integrated over 2/d,«,t- 

Our Bayesian analysis determines two parameters, the best fit values of the maximum 
and minimum deuterium abundances, yD.max an d 2/D,min = 2/D,max - w, compatible with the 
errors in the data and the assumption that any spread in the observed abundances, above 
and beyond that expecte d from the e r rors, is real and reflects the variable ISM depletion 



of deuterium onto dust. iHogan et al.l (119971 ) were mainly interested in the minimum 4 He 



abundance, Y min , and we are primarily interested in the maximum ISM D abundance, 2/D,max- 
The likelihood in the {^.max^} plane is maximized to find the best estimate of the maximum 
undepleted D abundance, providing a lower limit to the ISM D abundance (?/d,ism > 2/D,max)- 
We note that in the absence of observational errors, the depletion parameter, w, is restricted 
to < w < t/D,max, since |/D,max > 2/D,min- in Figures 2 - 4, w starts at zero and the dashed 
lines show the boundary, w = y^> 



Unax- 



2.3. Results 

As may be seen from Figured], the distributions of the observed D abundances for the LB 
and nLB LOS are very different. The LB yn values show little gas phase variation from LOS 
to LOS, in contrast to the nLB D/H ratios, whose variation accounts for most of the factor 
~ 4 spread in the observed ISM D abundances. Indeed, all 21 LB D abundances are consistent 
with no variation (within the observational errors) around a weighted mean of log(yp,LB) — 
0.19, corresponding to ?/d,lb = 1-5. Our Bayesian analysis for all five probability distributions 
is in agreement with this and yields w ~ and yr>, m ax ~ Z/D,min ~ 15. Our results for three 
distributions (top-hat, positive-bias, negative-bias), shown in Figure [2] illustrate this result. 
A lower bound to the deuterium astration factor may be inferred from the up per bound to 



the LB deuterium abundances, j/Dmax = Vdlb, and the estimate of y^p from iPettini et al. 
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(|2008|). For central values of log(y D ,LB) > 0.19 and \og(y DP ) = 0.45 f lPettini et all 120081 ). 
l°g(/D,LB) < 0.26, corresponding to /d,lb < 1.8 El, consistent with the successful Galactic 
chemical evolution models identified in SR T. This value is also consiste nt with a wide range 
of chemical evolution models discussed in IProdanovic &: Fields! (120081 ) with both low and 
high infall rates of nearly pristine gas, as well as with a variety of initial mass functions. 

In contrast to the LB, the nLB LOS do show large variations among the observed gas 
phase D abundances (see Figured]). It is thus expected that the Bayesian analysis of this 
data subset will find evidence for w > at a statistically significant level. This indeed is 
seen in Figures 0] (top-hat) and |5] (M-shaped bias) where our results are shown separately 
for the data of the LB (top panel), the nLB (middle panel) LOS, along with those for all 49 
LOS (bottom panel). For the 28 nLB LOS the best fit values, from a top-hat distribution, 
are t/D,max = 2.1 and w = 1.5. The top-hat Bayesian analysis for the complete FUSE data 
set (Figure HI bottom panel) finds maximum likelihood values of ?/D,max = 2.0 and w = 1.3, 
corresponding to ?/D,min = 0.7. The non-zero depletion parameter found for all 49 LOS is 
driven by the large variations in yn for the nLB LOS. This value of j/D,max corresponds 
to /d < 1-4, which, within the errors in the measurements of ?/dp and ?/D,max, is marginally 
consistent with the fiducial chemical evolution model discussed in SRT and is consistent with 



a range of models discussed in IProdanovic fc Fields! (120081 ). where even larger infall rates are 
required for some initial mass functions. On the other hand, when the complete FUSE 
data set is analyzed with the M-shaped prior suggested by the referee (Figure [5] bottom 
panel), the maximum likelihood is obtained for !/ Dimax = 1.8 and w = 1.0, corresponding 
to ?/D,min = 0.8 and /p < 1-6. This result is consistent with the fiducial model adopted 
in SR T, but not with some of the other Galactic chemical evolution models discussed by 
them (jSteigman et al.l 120071 ) . S imilar to the nLB ca s e, thi s result is also consistent with 
some of the models explored in IProdanovic fc Fields! (120081 ). where the observations allow 
for both high and low infall rates depending on the choice of the initial mass function. 
However, the more recent initial mass functions are only consistent with larger infall rates. 



Significant figures and round off: The Hi and Di column densities listed in Table 2 of iLinskv et al 



(|2006l ) are typically given to two significant figures (the integers in front of the decimal place simply count 
the powers of ten). These column densities are used to find log(?/D) = 5 + log N(D i) - log N(H i), where the 
errors in the logs of the column densities are added in quadrature. As a result, the values of log(t/D) (and 
their errors) are only known to two signific ant figures an d , so t oo are the values of uy> (and their errors) 
derived from them. However, in Table 3 of ILinskv et al.l (|2006f ) the values of y n are given to thre e (or 
more) significant figures. For our Bayesian analysis we use the data in Table 3 of ILinskv et al.l (|2006), but 
the results presented here are generally rounded to two significant figures. As an example, if we had first 
rounded the individual values of ?/d,lb and ?/dp to two significant figures, we would have found j/d.lb = 1-5 
and j/dp = 2.8, leading to fu < ?/Dp/yDXB = 1-9, in contrast to the slightly different value quoted here, 
which follows from log(fn) < log(?/Dp) - log(j/D,LB) = 0.26 and /d = 10 26 = 1.8. 
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Fig. 2.— Likelihood contours (68%, 95%, 99%) in the w - y D ,max plane for the 21 LB LOS for 
the top-hat (top panel), the positive-bias (middle panel), and the negative-bias probability 
distribution (bottom panel), where w = ?/D,max — £/D,min is the depletion parameter. The 
best fit values - for all three probability distributions - are at j/D imax = 1-5 and w = 
(2/D,min = 2/D,max), rounded off to two significant figures. The dashed line, w = ?/D,max, 
represents the bound above which the results, in the absence of errors, are unphysical since 
2/D,min = 2/D,ma X - w becomes negative. 

As anticipated from Figure [TJ the bottom panels of Figures H] and confirm that the 
variation among the gas phase D abundances for all 49 FUSE LOS is too large to be ac- 
counted for by the observational errors (to / at much greater than 99% confidence), for 
all adopted probability distributions. Indeed, e.g., for the top-hat probability distribution, 
the D abundances for all 49 LOS span a range of nearly a factor of three, from ?/D,max = 2.0, 
down to y D:min = y D ,max -w = 0.7. 

The effects of the choices of the priors on the likelihood distributions for ?/D,max and 
w are shown for the complete FUSE data set on Figure [3] and on the bottom panel of 
Figure [5j While the results for the top-hat and positive bias distributions are quite simi- 
lar, the negative-bias distribution is noticeably different. The best fit parameters ?/D,max> 
w and resulting astration factor /p are summarized in Table 1. Also presented in the 
table is, AlnJ2? max , the difference between the logarithms of the largest maximum likeli- 
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Fig. 3. — Likelihood contours (68%, 95%, 99%) in the w - yr>,msx plane for all 49 LOS for 
the top-hat (top panel), the positive- bias (middle panel), and the negative-bias probability 
distribution (bottom panel). The dashed line is w = 2/D,ma X ; see Figure [2j 

hood and a maximum likelihood of a given probability distribution. While the M-shaped 
probability distribution has the largest maximum likelihood value, so that Aln«if max = 
In j2f maXi M — In Jzf max ,bias, this distribution is closely followed by the top-hat and positive-bias 
priors, suggesting that 1.8 < ?/D,max < 2.0. The negative-bias prior, favoring large depletion, 
yields the poorest fit to the data. As may be seen from Figure [TJ the large errors in the 
data mask which of our priors might provide the best fit. Because the top-hat distribution 
requires the simplest assumption (no preference), and since its maximum likelihood is similar 
to that of the M-prior, we to adopt yo.max = 2.0 ± 0.1 as our best estimate for the maximum 
of the gas phase, ISM D abundances. 



3. Summary and Conclusions 



In the decades prior to the current era of precision cosmology th e primordial abundance 
of deu terium provided the only quantitative cosmological baryometer ( IBoesgaard fc Steigman 
19851 ). Although, at present, the deuterium abundance is only measured along seven, high- 
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Fig. 4.— Likelihood contours (68%, 95%, 99%) in the w - y D ^ plane for the 21 LB LOS 
(top panel), the 28 nLB LOS (middle panel), and all 49 LOS (bottom panel), using the 
top-hat probability distribution. The dashed line is w = yo.max (see Figure [2]). 



redshift, low-metallicity LOS to background quasars (jO'Meara et all 120061 : iPettini et al. 
20081 ). the inferred primordial D abundance, y^p = 2 .8 ± 0.2, is in excellent agreement 
with the non-BBN inferred baryon density parameter (jSteigmanl 120071 ) . In the post-BBN 
Universe deuterium is destroyed as gas is cycled through stars, so that comparing the abun- 
dance of deuterium in the ISM of the Galaxy with the primordial D abundance provides 
an estimate of the virgin fraction of the ISM (i.e., the amount of gas presently in the 
ISM which has never been cycled through stars), constraining models of Galactic chemical 



Table 1: Results for the different shapes of the adopted priors. 



Bias Shape 


1/D,max 


W 


/d 


A In .if 


Top-hat 


2.0 ±0.1 


1.3 ±0.2 


1.4 ±0.1 


0.5 


Positive 


1.9 ±0.1 


1.4 ±0.2 


1.5 ±0.1 


0.9 


Negative 


2.4 ±0.2 


1 7+0.3 
l - ' -0.2 


1.2 ±0.1 


2.7 


M-Shaped 


1.8 ±0.1 


1.0 ±0.1 


1.6 ±0.1 


0.0 


A-Shaped 


9 O+0.2 
z -°-0.1 


1 s+ - 3 


1.2 ±0.1 


1.4 
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Fig. 5.— Likelihood contours (68%, 95%, 99%) in the w - y D ,max plane for the 21 LB LOS 
(top panel), the 28 nLB LOS (middle panel), and all 49 LOS (bottom panel), using the 
M-shaped probability distribution. The dashed line is w = yo.max (see Figure 2). 



evolution (jSteigman fc Tosil Il992l . 1199.4 ISteigman et all 120071 : IProdanovic fc Fields! 120081 ). 
According to conventional wisdom the deuterium-free, metal-enhanced products of stellar 
nucleosynthesis should be well-mixed in the local ISM. In contrast, the FUSE data on the 
abundances of deuterium and several metals (e.g., iron, oxygen, etc.) along LOS within 



1 — 2 kpc of the Sun reveal a much different picture. The FUSE (ILinskv et al- 



and earlier observations flJenkins et 



al 



1999; Sonneborn et al. 2000: Hebrard et al 



20061 



2002: 



Hoopes et al. 2003 ; Prochaska et al. 20051 ) reveal unexpectedly large gas phase variations in 
?/d (and in the abundances of iron, oxygen, etc.) within the local ISM, as shown for FUSE 
deuterium data in Figure [TJ It has been proposed that the large variations observed in the 
local ISM D abundances can be accounted for by prefere ntial depletion of deuterium (relative 
to hydrogen) onto dust (1 Jural Il982l ; iDraind 12004 . 120061 ). although incompletely mixed infall 
of relatively unprocessed, deuteri um-enhanced, metal-fre e material may have contribu ted to 
some of the observed variations ( ISteigman &: Tosil Il992l . Il995t ISteigman et al.l 120071 ) . The 
large variations among the ISM D abundances, along with observational errors and the pos- 
sible contributions from dust depletion and infall, complicate using the D observations to 
provide a robust estimate of the ISM D abundance which, in combination with the primordial 
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D value, can lead to a constraint on the deuterium astration factor, /r> The key question 
is, given the data (with its errors), how to find the best estimate of the "true", undepleted, 
ISM D abundance? 

Here, to address this question, the limits to the true, undepleted, ISM D abundance 
were inv estigated employing a model-independent Bayesian statistical analysis similar to that 
used by iHogan et al.l (119971 ) to infer the primordial helium abundance f rom a set of helium 
abundance observations. It was assumed, along with lLinsky et al.l (120061 ). that the spread in 
the observed D abundances is the result of incompletely homogenized D depletion onto dust 
in the local ISM. In our analysis this is modeled by five different probability distributions 
(priors) for the yo values. The y^> (actually, log y^) values shown in Figure [1] suggest that, 
given the relatively large errors, a uniform (top-hat) distribution, favoring neither low-D nor 
high-D may be a good approximation to the data. To explore the sensitivity of our result 
to the choice of the prior, we first considered two asymmetric distributions - a positive-bias 
prior favoring low depletion, and a negative-bias prior favoring large depletion, as well as 
two other priors - an M-shaped distribution favoring both low and high depletion, and a 
complementary, A-shaped distribution. 



Using the FUSE deu terium observations along all 49 LOS (ILinsky et al.ll2006l ; lOliveira &: Hebrard 
20061 ; iDupuis et al.ll2009l ). we found the likelihoods in the {?/D,max,w} plane for the five choices 



of the Bayesian priors (see Figures [3] and [5]) . For all priors, the Bayesian analysis of the full 
data set requires significant depletion (e.g., w ^ at greater than 99.9% confidence). Com- 
paring the maximum likelihood values for the five different distributions, we find that the 
bimodal, M-shaped distribution provides the best fit to the observed data (see Table 1). 
However, it is important to notice that the shapes of the priors require an additional as- 
sumption in our analysis, so that the M-shaped distribution is the most model-dependent. 

1 

In contrast, the top-hat prior is the least model-dependent, favoring all levels of deple- 
tion equally. Given our ignorance of the detailed depletion mechanisms responsible for the 
observed scatter in the gas phase ISM deuterium abundances, we prefer to adopt for our 
estimate of the undepleted, ISM deuterium abundance, the result of the simplest, top-hat 



3 The M-shaped prior favors both low and high levels of D depletion while strongly disfavoring intermediate 
depletion, suggesting that two competing processes may be at work: depletion onto dust and evaporation 
from dust, perhaps due to exposure to shocks. To fit the M-prior scenario, both processes would have to 
be efficient and rapid to account for the deficit of intermediate D abundances. The distribution of the 
presently available data (with its errors) is inconclusive and does not strongly favor any of the adopted prior 
distributions. When more data become available, the Bayesian approach presented here may be used to 
learn more about the mechanism of deuterium depletion onto dust. 
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prior, whose maximum likelihood is similar to that of the best-fitting M-distribution, 



Zto.ISM > 2/D,max = 2.0 ± 0.1 = 2.0(1 ± 0.05) 



(10) 



This value is our best estimate of the true ISM D abundance based on the available deuterium 
observations in the local ISM and is independent of any model-dependent assumptions about 
galactic chemical e volution. Combining our result with y^p = 2.8 ± 0.2 = 2.8(1 ± 0.07) 
( iPettini et all 120081 ) (which, recall, provides a lower bound to the primordial abundance), 
yields a limit to the deuterium astration factor /b < 1.4 ± 0.1 (for the M-p rior, fn < 1.6 ± 



0.1), consistent with most, but not all, Galactic chemical evolution models (jSteigman et al. 



20071 ; iProdanovic &: Fields! 120081 ; iRomanol 120091 ). If, on the other hand we compared this 
IVrysM value to t he BBN + WMAP i nferred primordial D abundance, e.g., yp>p = 2.5 ± 0.1 
( ISteigmanl 120101 ). y^p = 2.5 ± 0.2 (jCyburt et all 120081 ) or, the pre diction inferred whe n 
including the WMAP-determined effective number of neutrino species ( IKomatsu et al.ll2010l ). 
Vdp = 3.0±0.4, the resulting deuterium astration factor would be somewhat lower in first two 
cases, fn ~ 1-3 ± 0.1, which is marginally problematic for some GCE models. In contrast, 
for the IKomatsu et all (120101 1 value of yr>p, fp> < 1.5 ± 0.2, which is entirely consistent with 
GCE models. 

As seen in Figures [H - HI for the LB there is little scatter among the gas-phase D 
abundances. The small scatter is entirely consistent with the observational errors (w = 0) 
and all LB D abundances are consistent, within the errors, with ?/d,lb = 1-5(1 ± 0.03). This 
suggests that for the Local Bubble, yp>,isM > 2/d,lb and, /d,lb < 1-8 ± 0.1, consistent with 
all the successful chemical evolution models identified in SRT. However, while the uniform 
LB D abundance suggests that D may be undepleted in the LB, for all LOS, ?/D,max ~ 1-3 
2/d,lb> suggesting either that D is depleted uniformly in the LB or, that outside of the LB the 
gas phase deuterium abundance may have been enhanced along some LOS by the addition 
of nearly primordial gas which has recently fallen into the disk of the Galaxy in the form of 
cloudlets which take some time to mix with the pre-existing gas in the ISM. Does ?/d,lb = 
1.5 or 2/D,max = 2.0 provide the best estimate of the lower bound to the ISM D abundance? 
If deuterium is depleted onto du st, why is there not a strong correlation between deuterium 
abundance and iron depletion^ (jLinsky et al. 2006 ) and, which refractory element is then 
best to use as proxy for determining deuterium depletion onto dust? These questions cannot 



4 As pointed out by the Referee, shock strength may play a role in accounting for the scatter observed 
in the correlation between the gas phase D and Fe abundances. If deuterium is loosely bound to the grain 
mantle while iron is locked into the core of the dust grain, deuterium would be more easily returned to the 
gas than iron when grains are exposed to shocks of modest strength, while iron might be removed from dust 
grains only by stronger shocks. The scatter in the correlation between the gas phase D and Fe abundances 
may be an indicator of shock strength. 
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be an swered by the analysis presented here. In a companion paper (jSteigman fc Prodanovic 



20101 ) . abundances of refractory elements are used in concert with the deuterium abundances 



in an attempt to resolve this question. 
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